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Abstract 

Living cells provide a fluctuating, out-of-equilibrium environment in which genes must coordinate 
cellular function. DNA looping, which is a common means of regulating transcription, is very much 
a stochastic process; the loops arise from the thermal motion of the DNA and other fluctuations 
of the cellular environment. We present single- molecule measurements of DNA loop formation and 
breakdown when an artificial fluctuating force, applied to mimic a fluctuating cellular environment, 
is imposed on the DNA. We show that loop formation is greatly enhanced in the presence of noise 
of only a fraction of k^T , yet find that hypothetical regulatory schemes that employ mechanical 
tension in the DNA-as a sensitive switch to control transcription-can be surprisingly robust due 
to a fortuitous cancellation of noise effects. 
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Genes operate within a crowded cellular environment that is constantly interacting with 
the encoding DNA through various proteins that bind along its contour. We are only be- 
ginning to understand how chromosomal packing, and a crowded and constantly fluctuating 
cellular interior affects genetic regulation |l| ; however, the cytoskeleton of eukaryotes already 
provides a startling example of the dynamic effects of the intracellular medium on cellular 
function [2|. The cytoskeleton is often considered an active polymeric gel, since its dynamic 
organization is driven by ATP hydrolysis, and is known to generate mechanical stress and 
shear forces. These active properties give rise to non-equilibrium effects such as the gener- 
ation of an "effective temperature" that can drive processes of embedded elements beyond 
the levels of thermal activation. Recently, force fluctuations an order of magnitude larger 
than thermal fluctuations have been measured within the cytoskeleton [3j. Such findings 
remind us that the cellular interior is very much an active, non-equilibrium medium, the 
effects of which need to be given careful consideration in the context of gene expression as 
well. 

A regulatory mechanism that is potentially acutely sensitive to environmental noise is 
the fluctuation driven formation of protein-mediated DNA loops. In genomic DNA, this is a 
ubiquitous motif for the transcriptional control of gene expression [4|. The lac operon, which 
is responsible for efficiently metabolizing lactose in E. coli bacteria, provides a canonical 
example of DNA looping. A lac repressor-mediated DNA loop is formed when tetrameric 
LacI protein simultaneously binds to two lac operator sites and is crucial for the repressive 
regulation of lac genes [sj]. Thermal fluctuations, which generate tiny entropic forces on the 
order of only fcgT '/l p pa 80 fN, where l p = 50 nm is the persistence length 6], are sufficient 
to form loops within the DNA, making the association rate of loops extremely sensitive to 
tension along the DNA molecule [?J ■ 

To explore the effects of environmental fluctuations on protein-mediated DNA loops, a 
1316 bp dsDNA molecule with two primary lac operators spaced 305 bp apart was tethered 
to a coverslip and then attached to an 800 nm polystyrene microsphere. The microsphere 
was then trapped within the linear region of the optical potential of a focused laser beam 
allowing us to apply a well defined tension to the DNA. Details of our axial-constant force 
optical tweezers as well as a discussion of the DNA preparation can be found in 8|. Ten- 
sion in the DNA was calibrated to include both the applied optical force, which is linearly 
proportional to the laser intensity modulated by an acousto-optic modulator (AOM), and 
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volume exclusion effects arising from entropic interactions between the microsphere and the 
coverslip 9j. The looped and unlooped states of the DNA molecule, which correspond to 
different axial positions of the microsphere, were measured by analyzing defocused images 
acquired on a CCD camera at 100 fps in the presence of 100 pM of LacI protein. This 
method provides excellent temporal resolution for detecting loop formation and breakdown 
events with time windows as short as 300 ms. 

Fluctuating forces were applied to the DNA by modulating the intensity of the trapping 
laser with an AOM connected to a data acquisition board and controlled by a custom 
Lab VIEW program. The program generated Gaussian white noise, at a sampling rate 
of 1/St, which was superimposed upon a set average optical force. The modulation was 
performed such that the force applied to the trapped microsphere was randomly chosen 
from a normalized Gaussian distribution of standard deviation a. 

It is the tension along the DNA that should display white noise statistics; however, by 
optically shaking the microsphere the resulting time correlation of the induced tension will 
effectively be low-pass filtered: 

fa(fMO> = -e-W r >, (1) 

where rj(t) is the DNA tension and a = a 2 5t is the noise strength. The characteristic 
time may be approximated by r c j/k, where 7 = 6nr]r is the hydrodynamic friction 
coefficent of the microsphere, r\ is the viscosity of the medium and r is the radius of the 
microsphere. If we only consider events that happen on timescales greater than r c we may 
approximate the colored spectrum of Eq. (JTJ) by purely white noise such that Eq. (JTJ reduces 
to lim Tc ^o(r](t)ri(t')) w 2a5(t — if). From the worm-like chain model for DNA, r c ranges 
from 5-8ms at tensions from 180-120 fN. Therefore, if we set 5t below the cutoff imposed 
by t c (for our experiments, we fix St at 2 ms (1/St = 500 Hz)) then the applied fluctuations 
are essentially white up to frequencies l/r c , and we may adjust the strength of the noise 
by tuning the width of the distribution a. It is more intuitive, perhaps, to consider the 
applied fluctuations as generating an additional energy contribution to the thermal modes 
of the DNA of magnitude ksT a = <j 2 5t/2 r j ) and to parametrize the fluctuations relative to 
the ambient temperature of thermal fluctuations T a /T. 

Figure [U shows the distribution of lifetimes for the looped and unlooped states at a mean 
tension of 153 fN and with fluctuations of T a /T = 0,0.01,0.05, and 0.12. As indicated 
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FIG. 1: Experimental measurements of the (a) looped and (b) unlooped cumulative probability 
distributions for various noise conditions (right to left: T a /T = 0, 0.01, 0.05, 0.12) at a mean tension 
of 153 fN. The solid lines are exponential (biexponential) fits to the looped (unlooped) data. The 
insert to (a) shows a typical raw trace of the DNA extension vs. time. 

in the figure, the lifetime of the looped state is independent of the fluctuations. This is 
consistent with previous findings that the looped state is insensitive to femtonewton forces 
\)\. The lifetime of the unlooped state, however, is clearly seen to decrease as we increase 
the fluctuations. At 2a, approximately 95% of the noise distribution is accounted for. Since 
volume exclusion forces are on the order of 35 fN [9j, a = 60 fN or T a /T = 0.12 are the 
largest fluctuations we can apply to the DNA without significantly clipping the distribution. 

We fit the cumulative probability distributions using the kinetic scheme detailed in [?| 
to extract loop dissociation and association lifetimes. In summary, the looped lifetimes are 
simply fit by a single exponential function parametrized by the looped lifetime r^: S(t) = l — 
exp(— t/ri). The unlooped kinetics, however, are more complicated, and may accurately 
be described by collecting all time intervals beginning with an unlooping event and ending 
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FIG. 2: Normalized unlooped lifetimes as a function of applied noise. The (square) data points 
were taken at a mean tension of h = 153 fN. The (diamond) and (circle) data points were taken 
at h = 123 and 175 fN respectively. Theoretical curves are shown for h = 123, 153, 175 fN. 

upon the formation of a loop: 

S 2 S S 1 L. (2) 

ki 

States Si and 5*2 arise because there are multiple unlooped sub-states available to the protein- 
DNA system. Si represents a state with only one occupied operator, which may either loop 
at rate to form state L, or remain unlooped and convert at rate ki to state S 2 . State 5*2 
is an alternate configuration with both or neither operator occupied, which cannot form a 
loop, but may convert back to state Si at a rate k 2 . The first-order kinetics results in the 
following biexponential function for the cumulative probability distribution: 

L(t) = 1 - ±- [c+e-^- - c_e-^+] , (3) 
la 

where k = k 2 + ki, a = [(k + k L ) 2 — 4:k 2 k L \ 1 ^ 2 , c± = (k — k L ± a) and the time constants are 
defined as r± = 2/(/t + k^ ± a). From this fit equation we are able to extract the unlooped 
lifetime r u = \jk^. In accord with the constant force results [7J], the interconversion rates 
k 2 and £q were found to be essentially independent of the applied tension. 

Figure [2] displays the unlooped lifetimes as a function of the noise T a /T normalized to 
the zero noise lifetime, T u (h;T a /T)/T u (h;0), about a mean applied tension of 153 fN. The 
results demonstrate that fluctuations do indeed drive loop formation and coalesce nicely 
with previous observations that femtonewton forces can radically affect the rate at which 
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Lacl-mediated DNA loops form [7|. Our results imply that an access contribution from 
environmental fluctuations of only 5% of the ambient thermal temperature T can double 
the rate at which DNA loops form, which is only a fraction of that which the cytoskeleton 
can induce upon an embedded polymer jjj. This stochastic mechanism might, therefore, 
provide an alternate 'noisy' means for mechanical control of genetic transcription. 

Although the rate at which DNA loops form is quite sensitive to environmental fluctua- 
tions, our data also show that this sensitivity is practically independent of the mean applied 
tension in the DNA. A separate measurement of the loop formation rate as a function of 
mean tension, h = 123, 153, and 175 fN, collected at a constant applied noise, T a /T = 0.05, 
reveals a striking observation: the normalized lifetimes, r u (h; 0.05) /r u (h;0), are constant 
with an approximate value of 0.5 irrespective of the average tension h in the molecule (see 
Fig. [2]). As we will show, this could allow regulatory schemes that are based upon protein- 
mediated DNA loops to display a significant level of robustness to noise. 

One might imagine that the increase in looping occurs because the fluctuations allow 
the random sampling of a low tension regime. However, any increase in looping gained by 
reducing the tension is negated by the equally likely sampling of high tension fluctuations, 
requiring us to consider the role of fluctuations in more detail. We begin with an effective 




FIG. 3: The loop association process can be modeled by diffusion over a barrier. The unlooped 
lifetime t u is given by the average time it takes for the DNA to diffuse from the equilibrium position 
xq, within the energy landscape, to the top of the energy barrier, of magnitude AU at x%, where 
it can escape to infinity and form a loop. 
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Langevin equation for the motion of the tethered particle system: 



dx dU (x; h) 



1 dt 



dx 



(4) 



The reaction variable x(t) describes diffusion along the energy landscape provided by U(x; h). 
The stochastic term £(t) accounts for thermal fluctuations and is modeled as a white noise 
source with zero mean, i.e., (£(£)£(£')) = 2A<5(£ — t'), with A = 7/csT. 
We choose the following phenomenological form for the energy barrier: 



U(x; h) = -ax 2 



1 „ 

-ox + hx. 
3 



(5) 



The harmonic term contains the cost of bending the DNA, while the cubic term, although 
somewhat arbitrary, is the simplest contribution that can give rise to an energy barrier 
between a second equilibrium state, which is here assumed to be the looped configuration, 
see Fig. [3j It should be noted that the parameter a is quite different from the spring constant 
k of a worm-like chain model since a must parameterize the full three dimensional search of 
the polymer to form a loop within this lower dimensional, effective energy landscape. The 
linear term represents the force h that we apply with optical tweezers to stretch the DNA, 
which effectively modulates the energy barrier by tilting the energy landscape and, therefore, 
increasing or decreasing the barrier height AU. Note that we are not attempting to account 
for the unlooping process, since the looped state, as our data reveal, is not sensitive to the 
forces we apply. 

An exact formula for the mean passage time 
is given by 

7 f 1 , f U(x;h) \ f x f -U(y,h) \ 

T " = YT ' 6XP l >- ^ I I 1 Z nn ) ■ ( 6 ) 



^ f U(x;h) 



10j across the energy barrier from Xq to X\ 
U(y;h)\ 



"B-L Jxo V ^B-L J V I^B-> 

If the potential barrier AU is large compared to /cgT, then Eq. (jH]) can be expanded about 
the vicinity of Xq and X\ to yield the Kramers formula for the unlooped lifetime 



27T7 



y /U"(x ;h)\U"(x 1 ;h)\ 
where AU = U(xi] h) — U(xo; h). From Eq. 

27T7 



exp 



AU' 



(7) 



exp 



we can rewrite this relation as 
(a 2 + Abhf' 2 ' 



6b 2 k B T 



(8) 



Va 2 + 4bh 

We must now determine the coefficients a and b that parameterize our model potential 
(Eq. ([5])). We do this by an iterative least-squares fit of the lifetimes, given by Eq. fl6]), 
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FIG. 4: Unlooped lifetime function of a constant applied force. The solid line is the 

theoretical fit, given by the Kramers relation of Eq. ([8]) to the data points (squares). The coefficients 
of the potential U(x; h) are a = 4.64 x 10~ 5 N/m and b = 739 N/m 2 . 

to our constant force (zero fluctuation) data, the result of which is shown in Fig. HI We 
have found the difference in the resulting fit parameters to be negligible between the exact, 
Eq. ([6]), and approximate Kramers relation, Eq. (jSJ). Moreover, the zero tension (h = 0) 
energy barrier AU = a 3 /Qb 2 = 7.4/c#T, from this fit, is in remarkable agreement with the 
energy cost for loop formation predicted by an elastic rod model of DNA, AU = 7.7fcsT 



Since we add noise to the system by linearly modulating the tension applied to the DNA, 
we may incorporate this additional noise by modifying the correlations of the stochastic 
source such that A = 7/cbTe, where we have introduced the effective temperature 



With the replacement T — > Te in Eq. (|Sj) we are able to account for the effects of noise 
on the unlooped lifetimes T u (h;T a /T) at a mean tension h. This model gives excellent 
agreement with our experimental measurements of the looping lifetime T u (h;T a /T) as a 
function of noise (see Fig. [2]). Furthermore, it predicts a robustness to noise similar to 
what we observe experimentally (see Fig. [2]). The theory predicts a fairly constant ratio 
r u (h; .05) / r u (h; 0) = 0.65,0.63, and 0.62 at a mean tension h = 123,153, and 175 fN, 
respectively. 
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T E = T(1 + TJT). 



(9) 
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Our results suggest how a force dependent genetic switch that employs DNA looping 
to regulate transcription could operate stably within a noisy environment. For instance, 
consider a regulatory element controlled by the formation of a DNA loop at a basal rate k tl 
under constant tension h\. A regulatory signal could be provided by a change in tension, h 2 , 
such that the loop formation rate is now a factor of p times the basal rate, i.e. k t2 /k tl = p. 
If we assume that the tension felt by the DNA fluctuates around the targets t\ and ti such 
that we have two new looping rates k tl and k t2 , our results imply that k tl /k tl = k t2 /k t2 , 
which means that k^/k^ = p, so that the expression signal is unaffected by the noise. 

We have demonstrated how noise from thermal and environmental fluctuations drives 
protein-mediated DNA loop formation, yet leaves the loops unaffected once formed. En- 
vironmental fluctuations only a fraction of the size of thermal fluctuations in the DNA 
can greatly enhance the rate at which these loops form. We interpret these results with 
a fluctuating barrier model that can quantitatively explain and predict our measurements. 
This model is based on the previously demonstrated sensitivity of loop formation to static 
mechanical tension, which led to the suggestion that cells may utilize tension to regulate 
transcription through mechanical pathways, as opposed to the more commonly considered 
biochemical ones [7|. Based on our new observations, we may now postulate the feasibility 
of an alternate mechanical regulatory mechanism that uses environmental fluctuations as a 
means to control transcription. The rate enhancement of fluctuations might, in fact, explain 
why a simple loop results in several hundred fold repression in vivo 5|, even though the 
looped and unlooped states have roughly equal lifetimes in vitro. Furthermore, we have 
shown that the sensitivity of loop formation to fluctuations is insensitive to baseline static 
mechanical tension, and have demonstrated how this feature can lead to a robustness in 
regulatory function. 
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